Spin-charge separation in two-component Bose-gases 
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We show that one of the key characteristics of interacting one-dimensional electronic quantum 
systems, the separation of spin and charge, can be observed in a two-component system of bosonic 
ultracold atoms even close to a competing phase separation regime. To this purpose we determine 
the real-time evolution of a single particle excitation and the single-particle spectral function us- 
ing density-matrix renormalization group techniques. Due to efficient bosonic cooling and good 
tunability this setup exhibits very good conditions for observing this strong correlation effect. In 
anticipation of experimental realizations we calculate the velocities for spin and charge perturbations 
for a wide range of parameters. 



One of the most exciting recent events is the ever- 
growing interplay between previously disconnected fields 
of physics, such as quantum optics and condensed matter 
physics. In particular, cold atomic systems have opened 
the way to engineer strongly interacting quantum many- 
body systems of unique purity. The unprecedented con- 
trol over interaction strength and dimensionality allows 
the realization of "quantum simulators" where funda- 
mental but hard to analyze phenomena in strongly corre- 
lated systems could be observed and controlled. Exam- 
ples are the observation of superfluid to Mott insulator 
transition for Bose gases [ij and the fermionization of 
strongly interacting one dimensional bosons \% [3[ . 

Among interacting systems, the physics depends very 
strongly on the dimensionality. In one-dimensional sys- 
tems the interactions play a major role and lead to dras- 
tically different physics than for their higher dimensional 
counterpart. Typically, interactions in one-dimensional 
systems lead to a Luttinger liquid state where the exci- 
tations of the system are collective excitations The 
importance of such a state for a large variety of exper- 
imental devices in condensed matter has led to a hunt 
to observe its properties. A remarkable consequence of 
such a state is the absence of single particle excitations. 
This means that a quantum particle, that would normally 
carry both charge and spin degrees of freedom, fraction- 
alizes into two different collective excitations, a spin and 
a charge excitation. Such a fractionalization of a sin- 
gle particle excitation is the hallmark of collective effects 
caused by interactions. However, just as detecting frac- 
tional excitations in the case of the quantum hall effect is 
difficult , observing spin-charge separation has proven 
elusive despite several experimental attempts [1, 0, 
So far, the best experimental evidence is provided by 
tunneling between quantum wires where interferences ef- 
fects are due to the existence of two different velocities 
0. However, in these systems it is hard to quantify or to 
tune the interaction between the particles which causes 
the collective effects. Since control of the interaction is 
a routing procedure in ultracold gases, the possible real- 
ization of the phenomenon of spin charge separation has 



also been discussed in the context of cold fermionic gases 
IC, U, [13, and strongly interacting bosonic gases 

a. 



However, proposals to observe spin-charge separation 
in ultracold fermionic gases are still plagued by the cur- 
rently quite high temperatures in such systems. A much 
better setup to test for spin-charge separation would 
be two-component Bose gases, for example using the 
\F — 2, nip = —1) and the \F = l,mF — 1) hyper- 
fine states of ^^Rb [l^ [l^. Experimentally, this sys- 
tem retains the advantages of the fermionic ultracold 
atom setup while allowing for much lower temperatures 
due to the more efficient cooling techniques available for 
bosons. However, theoretical studies [13, [11] for one- 
dimensional systems predict that close to the experimen- 
tally accessible parameter regime of almost equal inter- 
and intra-species interaction strength phase separation 
occurs. This is the remaining potential experimental 
complication in the setup. 

In this work we demonstrate the phenomenon of spin- 
charge separation in the experimentally relevant parame- 
ter regime, allowing to use this system to unambiguously 
test for spin-charge separation. We calculate both the 
real-time evolution of a single particle excitation and the 
dynamical single particle spectral function of the two- 
component bosonic systems. We show that both these 
quantities demonstrate the separation of a single-particle 
excitation into spin and charge. We further determine 
the velocity of spin and charge and the Luttinger param- 
eters for experimentally relevant parameter regimes. To 
perform the calculations we use variants of the density 
matrix renormalilzation group method (DMRG) |1S(,, JU] ■ 
The numerical treatment is necessary to obtain reliable 
predictions for experimental realizations, due to the close 
proximity of this regime to phase separation. 

A one-dimensional two-component Bose gas in an op- 
tical lattice filll can be described by the two-component 
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Bose-Hubbard model 



(1) 

Here j is the site index and 1^ — 1,2 labels the two dif- 
ferent hyperfine states of the system, b and 6^ are the 
annihilation and creation operators and fi is the num- 
ber operator. The first term models the kinetic energy 
of the atoms. The intra-species interaction is described 
by the Ui, term. We use U := Ui = U2 as it is ap- 
proximately realized for commonly used hyperfine states 



of ^ Rb [22|. The inter-species interaction is given by 



2 and the 'spin' density 
- nj,2- We focus on systems with average fill- 
rij^i/L = nj^2/L smaller than one particle 



the U12 term and the last term describes external po- 
tentials. In the following we use the dimensionless pa- 
rameters u = U/J and U12 = Uu/J- We define the 
'charge' density n-,- c = Uj^i + n 

ing n = J2j 

per site and parameter regimes up to close to the tran- 
sition to the phase separation (approximately U12 ~ u 
[3]). Here L is the number of sites in the system. 

In a superfluid phase away from the transition to phase 
separation the low energy physics can be approximated 
by a density-phase representation of the bosons as used 
in the bosonization method [j]. In this representation 
the Hamiltonian is totally separated into one part for 
the charge and one for the spin degrees of freedom. The 
physics is fully determined by the velocities Vc^s and 
the so-called Luttinger parameters Kc^s for spin (s) and 
charge (c). Therefore the separation of a single particle 
excitation into spin and charge excitations is expected. 
The parameters of two interacting species of bosons can 
be related to the parameters K and vq for the single 
species case [l^l by 



Vc,s = vo\/l ± {gi2K) / (ttvo) 
and Kc,s = K/^l±igi2K)/iTTVo). 



(2) 



Here gi2 is the interspecies interaction strength in the 
continous model. In the limit of small interactions, the 
single species parameters K and vq can be directly re- 
lated to the Bose-Hubbard Hamiltonian Eq. ^ f3]. For 
higher values of the interaction strength the relation even 
for the single species situation is unkown, and has to 
be determined numerically [24] . For large values of the 
inter-species interaction of the order of the intra-species 
interaction the system approaches the transition to phase 
separation and the bosonization approach becomes a pri- 
ori inaccurate. 

Snapshots of the real time-evolution of a single par- 
ticle excitation in a two-component bosonic system are 
shown in Fig. [TJ The single particle excitation at time 
t = is prepared by the application of the creation op- 
erator of one species, say 1, on site L/2 to the ground 
state \tpo), i.e. b^j^^^ ilV'o)- The resulting sharp peaks in 



the density distributions are shown in Fig. [T] (a). For 
t > the time-evolution of the single particle exci- 
tation is calculated using the adaptive time-dependent 



DMRG [25|, |2a] ■ The time-evolution is performed using 
a Krylov algorithm 27| in a matrix product state basis 
with a fixed error bound for each timestep of the or- 
der of IW-^it + At)) - e:iq)[-iHAt]\'^{t))f < IQ-^. The 
stepsize At — 0.2 and 6 to 10 Krylov vectors were used 
resulting in Hilbert spaces with a local dimension of a 
few thousand states. As can be seen in the snapshots in 
Fig. [T] the initial single particle excitation splits up into 
two counter propagating density waves. Due to their dif- 
ferent spin and charge velocities, after a period of time 
a clear separation of the density waves is seen (cf. Fig. [T] 

(c)) m- 

Additionally to the time-evolution of a single parti- 
cle excitation, we obtained the single particle spectral 
function A^{q,w) = ^^{b^^ fj^^^^^_j^^ bq^u) as shown 
in Fig. [21 For fermions this function is known to ex- 
hibit two pea ks at the spin and charge excitation en- 
ergies [2^ [2^ , showing a direct signature of the spin- 
charge separation. For the bosons computing this spec- 
tral function is more involved and up to very recently 
it was only derived for a single component bosonic sys- 
tem Isil . [3^ . An expression for the correlation func- 
tions which allow to obtain the single-particle correlation 
function within the bosonization treatment for a two- 
component bosonic system was derived in (ssj . Power 
law singularities at qVc.s are obtained with respective ex- 
ponents + l/2i^s,c — 1. For the values of the 
Luttinger parameters (as shown in Fig. [2]) one thus ex- 
pects two divergent peaks. We show in Fig. [2] the full 
spectral function for our microscopic model, as calculated 
numerically using a matrix product state generalization 
of the correction vector method 3l|, . Our results show 
clearly the appearance of the two separated peaks, the 
lower representing the spin and the upper the charge ex- 
citation branch 4^] . Thus both the real time evolution of 
a single particle function and the single particle spectral 
function show clear signatures of separation of spin and 
charge. 

To observe the separation of spin and charge exci- 
tations experimentally in a system of ultracold bosons, 
knowledge of the spin and charge velocities is indispens- 
able. We therefore determined the velocities for a wide 
range of parameters. This was done calculating the time- 
evolution of a small spin and charge density perturbation 
respectively. The density perturbation was created at 
time t = applying an external potential of Gaussian 

form e,y,j = eo^xp ^— ^^^^"^ ^ where v = c,s. At time 
t — Q the potential is switched off and the time-evolution 
of the density perturbation is calculated. The errors in 
the obtained velocities are of the order of 0.01aJ/?i for 
small U12 and increase with larger U12. 

In Fig. [3] we show the dependence of the velocities 
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FIG. 1: (Color online) Snapshots of the time-evolution of the charge and spin density distribution of a single particle excitation 
created at time t = Oh/J; (a) at time t = Oh/J, (b) at time t — 1.5h/J and (c) at time t — 2.5R/J. The system parameters 
were ni,2 = 0.625, u — 3, ui2 — 2.1. The charge density is shifted by 1.25 for better visibility. The arrows in (c) mark the clear 
separation of the charge and the spin density waves. 




FIG. 2: One-particle spectral function at momentum q = 
20/657r/a. Two peaks corresponding to the spin and the 
charge excitation can be distinguished. The following param- 
eters were used n = 0.625, u = 3, ui2 — 2.1 on a system with 
L = 64 sites and a broadening rj = 0.1. Here a is the lattice 
spacing. 



on the inter-species interaction for two different values 
of the parameter 7 — u/{2nc). For both parameter 
regimes, the charge velocity increases with increasing in- 
teraction whereas the spin velocity descreases. For a 
vanishing inter-species interaction it was shown in 
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that for 7 < 1 the bosonization and the solution of the 
exactly solvable Lieb-Liniger model agree approximately 
with the velocities in the Bose- Hubbard model. In Fig. [3] 
(a) we find for 7 « 1.1 very good agreement with the 
analytical solutions even up to strong inter-species in- 
teraction strength U12 close to phase separation. For 
7 « 2.4 [cf. Fig. [3] (b)] even for vanishing interspecies in- 
teraction the deviations from the direct relation between 
the parameters in the Bose-Hubbard model are consid- 
erable (cf. [Si)- However, we find that the dependence 
of the velocities on the interspecies interaction strength 
via Eq. ( [2]) is still a good approximation correcting by 
the numerically determined value for vq and K. This 
holds even up to close to the regime of phase separa- 
tion, i.e. U12 ~ u where the difference in the velocities is 
maximal. However, for U12 « u the results for the spin 
velocity start to deviate for both values of 7. 



In Fig. [?] we show for two different fixed inter-species 
interaction strengths the dependence of the velocities on 
the density. The charge and the spin velocities rise with 
increasing background charge density. (Note, even at 
ric — 1 the system is in the superfluid regime.) The 
increase of the velocities is described to good accuracy 
using the analytical form Eq. provided we use nu- 
merically obtained values of K and vq- For large U12 
and small ric the results for the velocities from DRMG, 
in particular the spin velocity, deviate considerably from 
Eq. (12), showing that the approximate relation cannot 
be used in this regime. Note that in this regime the 
extraction of the spin velocity from the real-time evolu- 
tion becomes also more involved since the spin pertur- 
bation shows a strong spreading in space (cf. [35|). At 
the timescales over which we calculated the velocity, the 
left- and the right-moving spin perturbations are not yet 
fully separated and show strong amplitude damping. Our 
finding of the dependencies can be used to predict the ve- 
locities for experimentally interesting parameter regimes. 

In recent experiments for the preparation of a mixture 
of two bosonic components in optical lattices mostly two 
hyperfine states of ®''Rb are used, e.g. the \F = 2, nip = 
— 1) and the \F = l,mp = 1) hyperfine states. The 
intra-species scattering lengths are 02 = 91.28aB and 
ai = 100. 4aB (2^ . respectively where as is the Bohr ra- 
dius. For these states the inter-species scattering length 
is of the same order of magnitude as the intra-species 
scattering length and can be tuned about 20% using a 
Feshbach resonance 15, [l3|- Thereby the experimen- 
tal parameters are close to the competing phase sepa- 
ration regime. These mixtures can be confined to one- 
dimensional structures using strongly anisotropic lattices 
BBS, 37 1 . The most intuitive observation of the phe- 
nomenon of spin-charge separation in these systems is 
to generate a single particle excitation and then follow 
the evolution of the excitation in real time. This can be 
done measuring the spin-resolved density over a certain 
region. The creation of a single particle excitation can 
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FIG. 3: (Color online) Dependence of the charge and spin 
velocity on the interparticle interaction strength for (a) un = 
2 and n « 0.88 and (b) «i2 = 3 and n — 0.63. A comparison of 
analytical results (line, see text) and numerical DMRG results 
(symbol) is shown. The velocities are measured in units aj /%. 
Note that the errors of the DMRG results increase close to 

U-12 ~ U. 
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FIG. 4: (Color online) Dependence of the charge and spin 
velocity on the charge background density., comparison ana- 
lytical results of bosonization and numerical DMRG results. 
The parameters used are (a) u = 3, ui2 = 1.2 and (b) tt = 3, 
ui2 = 2.7. 



be done e.g. using outcoupling of single particles by the 
application of a magnetic field gradien t for addressability 
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The efficiency of 



and a microwave field 
such a technique for generating single particle exciations 
was demonstrated |39| using a cavity. The microwave 
field could be chosen to couple the \F = l,mp = 1) hy- 



perfine state to e.g. the \F = 2,mp = 2). This has the 
advantage that scattering with the |F = 1, = 1) state 
are supressed. The measurement of the density resolved 
over a region of approximately 10 lattice sites can then be 
performed using again the magnetic field gradient to get 
an unambigious signal. In an array of one-dimensional 
tubes the broadening of the signal caused by the trap- 
ping potential could be suppressed by preparing most of 
the tubes in a Mott-insulating state as shown in [l^ . 
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